Analysis of an atom laser based on the spatial control of the scattering length. 
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In this paper we analyze atom lasers based on the spatial modulation of the scattering length 
of a Bose-Einstein Condensate. We demonstrate, through numerical simulations and approximate 
analytical methods, the controllable emission of matter-wave bursts and study the dependence of 
the process on the spatial shape of the scattering length along the axis of emission. We also study 
the role of an additional modulation of the scattering length in time. 
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I. INTRODUCTION 

Atom lasers are sources of coherent matter waves that 
use an ultra-cold gas of trapped alkalii atoms as a reser- 
voir from which coherent pulses of thousands of atoms 
are extracted. Since the experimental achievement of 
Bose-Einstein Condensation (BEC) in gases pj, several 
methods have been proposed to deliver atoms from their 
confinement. The first of these devices used short radio- 
frequency pulses as outcoupling mechanism, flipping the 
spins of some of the atoms to release them from the 
trap [2| . Later, other atom lasers were built leading to 
ulsed, semi-continuous or single-atom coherent sources 



The quality of an atom laser is given by the amount 
of atoms that can be delivered from the trap and by the 
purity of the emission process. Concerning this point, it 
has been recently shown |8j that spin-flipping techniques 
present serious limitations in the number of atoms that 
can be emitted. This effect is related with the fluctua- 
tions at high flux due to the fact that the output-coupling 
mechanism populates all accessible Zeeman states. This 
constitutes a strong drawback for practical applications 
of atom lasers in high precision measurements like matter 
wave gyroscopes 

On the other hand, an atomic soliton laser using the 
mechanism of modulational instability was proposed in 
Ref. ^(j. I n this case, the emission is obtained by the 
combination of purely nonlinear effects in the atom cloud 
and the relaxation of the trap, it being necessary that the 
total number of particles in the cloud exceeds a critical 
threshold. This system has the advantage over standard 
atom lasers of pro ducing matter wave pulses in the form 
of solitons [ll|, Il2l Il3| , a kind of nonlinear waves gener- 
ated by a perfect balance between dispersive and nonlin- 
ear effects, yielding to robust wave packets that propa- 
gate without shape distortion [TJ| . Other types of atom 
laser based on nonlinear effects have been proposed [l5| . 

However, although one could extract a few coherent 
solitons from a Bose-Einstein Condensate by the mecha- 



nism of modulational instability, the final output would 
be very limited since the number of atoms per pulse gen- 
erated by this phenomenon is only a small fraction of the 
initial number of atoms in the condensate due to collapse 
processes. Moreover, the number of solitons generated is 
not large and half of them would be directed backwards. 
Finally, the trap must be destroyed for outcoupling and 
the pulses travel at different speeds once the trap is re- 
moved. Thus, it is important to discuss new outcoupling 
mechanisms for atom lasers. This is specially interesting 
since the techniques for generating BECs with large num- 
bers of particles and their physical properties are nowa- 
days well established and the current challenges in the 
field face the design of practical devices 0] . 

In a recent work 17 a novel outcoupling mechanism 
for an atom laser has been proposed. The method is 
based on the fact that a spatial variation of the scat- 
tering length (a) (see Fig. can be used to extract a 
controllable train of up to several hundreds of atomic soli- 
tons from a BEC without altering the trap properties. In 
Ref. 01 > a simple model for the spatial variation of the 
scattering length (a step function) was used in order to 
introduce the basic ideas and illustrate the phenomenon. 
However, for the practical realization of such device, a 
deeper analysis must be done, including studying the role 
of more physical distributions of the spatial dependence 
of a, and the temporal control of the output by the con- 
sidering an additional time dependence of the scattering 
length. In this paper we consider these aspects of the 
soliton emission process which are relevant for the the- 
oretical understanding and experimental demonstration 
of this new type of atom laser. 

The structure of this paper is as follows: in Section ITU 
we present the configuration of the system and the math- 
ematical model to be used for the theoretical analysis. In 
Section IlIII we study the mechanism of emission of single 
solitons. To do so we first study a simplified model for 
which an analytical study can be done. We also analyze 
theoretically a more realistic scheme by means of varia- 
tional methods and numerical simulations, for Gaussian- 
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shaped spatial variations of the scattering length. Fi- 
nally, in Section HVI we study the effect of temporal vari- 
ations of the scattering length induced by pulsed beams. 



II. SYSTEM CONFIGURATION AND 
THEORETICAL MODEL 

As in Ref. 0]> we assume a BEC which is strongly 
confined in the transverse directions (x, y) and weakly 
along the longitudinal one (z) leading to a cigar-shaped 
configuration. We will consider the effect of a spatial vari- 
ation of the scattering length along z from positive (or 
zero) to negative values. In principle this can be done by 
magnetic [fa | or optical fl9j means, for instance by using 
an appropriate laser beam that shines one of the edges of 
the condensate as in the sketch of the system plotted in 
FigH In this case the region of negative scattering length 
can be varied and displaced along the condensate by sim- 
ply moving the laser beam. The variation of a can also be 
controlled in time by using pulsed laser beams. Similar 
arguments apply for magnetically controlled scattering 
length, although the control of laser allows for faster and 
easier manipulation of the spatial variations of the scat- 
tering length. In any case what it is important is that 
we will depart from the model of step-like spatial vari- 
ation of the scattering of Ref. 0] to more realistically 
achievable smoother dependences. 

When the zone in which the scattering length is man- 
aged to negative values overlaps with the wings of the 
atom cloud, as in Fig^ the BEC feels a repulsive force 
whose strength depends on the number of atoms in the 
condensate. When this effect is large enough, the trap 
is overcome and part of the cloud is delivered and emit- 
ted outwards. When the condensate refills the gap left 
out by the outgoing pulse the process starts again and 
a new soliton is emitted. This process would continue 
while there is a large enough remnant of atoms in the 
trap and would lead to a soliton burst escaping from the 
BEC. 

Therefore, we will consider that the cloud of N equal 
bosons of mass m is tightly trapped in (x,y) by a har- 
monic potential V± of frequency v±_ and weakly confined 
along z by the effect of an optical dipole trap 14(20, 21] 
that can be produced by a laser beam of a given width 
along z (see Fig. Thus, we have 




FIG. 1: [Color online] Sketch of the system we will study in 
this paper for the case of optically controlled scattering length 
showing the BEC in the optical dipole trap, the transverse 
magnetic confinement and the laser beam used to manage 
the scattering length region. For a critical number of atoms 
or equivalently below a critical value of the scattering length, 
a burst of matter-wave solitons is emitted along the weakly 
confining axis. 



by a Gross-Pitaevskii equation of the form: 
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V 2 * + V(r)^> + U(z)\y\ 2 ^, (2) 



where $ is the condensate wavefunction, and its norm 
N = J 1^1 2 d 3 r gives the number of particles. The co- 
efficient U(z) = 4irh 2 a(z)/m depends on the scattering 
length a, which characterizes the 2-body interactions be- 
tween atoms. As commented previously, we will consider 
a to be a function of z with a localized region in which 
it becomes negative. 

We will concentrate on situations in which the spa- 
tial size of the ground state of the optical dipole trap 
is much larger than the ground state of the transverse 
harmonic potential leading to effectively one-dimensional 
dynamics. In this situation we can describe the dynam- 
ics of the condensate in the quasi-one dimensional limit 
as given by a factorized wavefunction of the form pl| 
*(r,t) = $o(x,y) ■ ip{z,t), satisfying 



dtp 



d 2 ^ 



2 dz 



(3) 



V(r) =V X + V, 
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f-exp [ - — 



(1) 



where Vb is the depth of the shallow optical dipole po- 
tential and L its characteristic width along z. The choice 
of a shallow Gaussian trap is very important, since the 
goal is to obtain outcoupling of solitons along the z axis. 
This can be achieved with a potential barrier that can be 
overcome by the self-interaction effects. The dynamics of 
the previous system in the mean field limit is described 



where rj_ = y/ h/mv± is the transverse size of the cloud, 
f(r)) — V z /(hi>j_), t — is±t is the time measured in units 
of the inverse of the radial trapping frequency and g(z) = 
2nr 2 L a(z) is the effective interaction coefficient. 

As an example we will present specific numbers in this 
paper corresponding to 7 Li, using the experimental pa- 
rameters of Ref. V = Hvx/2, v±_ = IK Hz, L = 4rx, 
N = 3 • 10 5 , w = 5Ar± a = —lAnm and times ranging 
from t = to t = Is. However, our results hold for dif- 
ferent atomic species like 85 Rb and 133 Cs, with adequate 
parameters. 
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FIG. 2: [Color online] Dependence of the number of particles 
in the condensate N on the chemical potential n for the simple 
case of square-shaped potential. The values are normalized 
by N m in and fMnin which correspond to the minimum of the 
curve. For the experimental values considered in the text for 
7 Li we find N m in = 2 x 10 4 and /imin = hu±/2 for a — —0.3 
nm. The dotted part of the curve (i(N) corresponds to an 
unstable region. The insets show the eigenmode profiles for 
different values of \i, the dotted lines indicating the bound- 
aries of the trap. The region of negative scattering length is 
the shaded one. 



III. SINGLE SOLITON EMISSION 



A. Exact theory for square potentials 



73 = 2fj,/r±hv±, we get 

{A 1 exp(-7iz) 
A 2 cos[(z - z 2 ) 72 ] 
A 3 sech[(z - 2:3)73 



z < -L 
-L < z < L 
z > L 



(5) 



As it can be appreciated in the insets of Fig. [3 the shape 
of the wavefunction depends on a*. Thus, a potential 
connected to a region where atom-atom interactions are 
activated, displays a continuum of stationary fundamen- 
tal states. For values of a* close to zero, (or equivalently 
a low number of atoms N in the BEC), the cloud is lo- 
cated at the center of the trap. As the product Na* is 
increased, the center of the cloud is displaced towards the 
boundary between zones two and three. For higher values 
of iVa 4 , the cloud is completely located in the region with 
a < and takes the form of a soliton with a hyperbolic 
secant profile. The continuity of the wavefunction of the 
condensate and its derivatives at the boundaries between 
the three regions yields to the following relationships: 



Z 2 = 

Z3 = 

A x = 

A 2 = 

A, = 



L H tan 

72 

L tanh 1 

73 



72 

73 



-71 

72 

tan [72 (L - z 2 )] 



A 3 sech[ 73 (L - z 3 )] sec [ 72 (X - z 2 )] 
A t cos[- 72 (i + z 2 )] e^ L 



sechz \ 



2sech 2 z - + 1 



7rrj_a* 



(6a) 

(6b) 

(6c) 
(6d) 

(6e) 



As we have pointed before, we will consider a situation 
in which the scattering length is changed on a localized 
region near the edge of the optical trap. In that case, it is 
possible to understand the basics of the phenomenon of 
the emission of a single soliton using a simple model. We 
consider a square trapping potential of width L and depth 
Vq and a step-like scattering length, as it is illustrated in 
Fig. El Mathematically 



V(z) 



z < —L 

Vq < -L < z < L 

z > L 



(4a) 



and 



a{z) 







z < L 



a* < z > L 



(4b) 



Although this model is a rough simplification of the more 
realistic situations to be described in detail later on, it 
captures the big picture of the process. In fact, this model 
can be solved analytically, by simply calculating the so- 
lution of the Schrodinger equation in the three zones. 
Defining 71 = 2fi/r±hi i ±, j 2 = 2(/i — Vo)/r±Hv± and 



The number of atoms N in the stationary state can be 
easily calculated by integrating \ijj\ 2 over z, 



A 2 

N=- 2 - exp( 
271 



A 2 

-2 7l L) + {1 



tanh [73 (L - z 3 )}} 



Aj 
2 



2L 



73 

sin (2j 2 L) cos (27 2 2 2 ) 



72 



(7) 



The previous equation provides the dependence of the 
chemical potential of the ground state on the number of 
atoms. As it can be seen in Fig. [2Jthe curve N vs /j, for all 
the stationary ground states has one maximum and one 
minimum, and displays a negative slope between them. 
It is well known from the theory of nonlinear Schrodinger 
equations [22ll23| that the stationary states correspond- 
ing to the zone with dN/dfj, < are unstable, what is 
known as the Vakhitov-Kolokolov's criterium of stabil- 
ity. This fact implies that any small perturbation will 
affect dramatically the eigenstates in this zone, that are 
basically nonlinear surface waves located in the bound- 
ary between the linear potential and the nonlinear region. 
As a consequence of this instability, the surface wave is 
reflected by the boundary and emitted to the nonlinear 
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B. Approximate theory for realistic potentials 

The analysis of the previous section, although exact, is 
only a crude representation of the real potentials which 
can be used in experimental scenarios. To study more 
physical spatial distributions of a we will combine the 
use of numerical methods with approximate averaged La- 
grangian formalisms p3 | consisting on minimizing the 
Lagrangian density over a family of trial wavefunctions 
which we will choose as 



ip(z, t) — Aexp 



(z-zo(t)Y 
2w 2 



exp[iv(r)z], (8) 



This ansatz leaves a single free degree of freedom for the 
atom cloud: its center zq(t), which moves with speed 
proportional to v(t). 

To study an scenario close to a possible experimental 
setup in which the scattering length is optically managed, 
we will analyze the soliton emission when the laser beam 
used to manage the scattering length takes the form 



a(z) 



aexp 



z-D 



z<D + w g (9) 

Z > D + Wg, 



i. e. an hypergaussian ramp connected at its maximum 
to an infinite plateau. This choice allows us to cover 
many cases between the two limits m = 2 corresponding 
to a Gaussian distribution and m — > oo representing a 



step function. The standard calculations taking a(z) as 
given by Eq. @ lead to a Newton-type equation of the 
form 



dzo 



(10a) 



where TL(zq) is given by 



II(zo) = r 2 



+ 



N 

2w 2 TT 



1 - 



exp 



1 + 1 
-2(z-z ) 



a{z)dz 



(10b) 



From Eqs. llOal) we see that the center of gravity of 
the cloud (zo) behaves like a classical particle under the 
effect of a potential n(z ). This provides a qualitative 
understanding of the soliton emission: for the linear case 
(a = 0) the center of the cloud is located at the bottom 
of the Gaussian trap and approximates its fundamental 
eigenstate [Fig. El (a)]; as a takes more negative values 
the effective trapping of the cloud is deformed and the 
minimum of the equivalent potential moves to the region 
with z > 0. 

The calculations, in the particular cases of m = 2 -a 
Gaussian beam- and m — > oo -a step function-, lead us 
to the effective potentials: 



n 9 (* ) = A -Ml- 




1 -I- ^- 

1 "I- L 2 



aN 
w\J2ti 



erfc 



V2(D- 



zq ) 



=exp 



{D-z y 



ZzpWg + 2w g D + 2w 2 g + w 2 \ _ ^ f ( -2z a w 2 g - Dw 2 



w^2w 2 + w 2 



w g wJ2w 2 + w 2 



Vn 



n s( , ) = V»i<^|i- 



1 aN 
/2n w 



erfc 



V2(D 



zo) 



w 




(11a) 
(lib) 



where g and s refer respectively to the step and Gaussian ramps. The functions erfc(u) and erf(w) = 
Jq exp(— v 2 )dv = 1 — erfc are the complementary error and the error functions, respectively. Fig. [3] shows 

the equivalent potentials II given by Eq. I|lla|) and Eq. lllbl) for different values of a. As the scattering length 
becomes more negative, it reaches a limiting value a cr for which the potential 11(0) = II(oo), thus if the atom cloud 
is initially placed at z = it will oscillate around the minimum and escape z (t) — > oo for r — > oo, a phenomenon 
which is called soliton emission [2!| . The critical value of a that corresponds to the threshold for soliton emission can 
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be obtained within our formalism from the condition, 11(0) = II(oo) which leads to: 



Nat 
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hf±Ji 



-^erfcf D + W °V2 



, 2w„D + 2wi + w 2 
erf | — 

Wt/2Wg + w 2 



erf 



ftv±_ \/L 2 + w 2 



erfc[^}-2 



2w 2 + w 2 



= exp 



-Dw 



w g j2w 2 + w 2 



-2D 2 



2w 2 + w 2 



2 

V2 



(12a) 
(12b) 



In both cases, it is possible to find approximate expres- 
sions that hold for D w w and Vq = Ahv±, with A < 1, 
which provides a shallow trap. Thus, it is possible to 
write: 



Nat 



-2.655 



-2.565 



AwL 



VT 2 



AwL 



\JL 2 + w 2 ' 



(13a) 
(13b) 



We have found a good agreement between the exact nu- 
merical value and the prediction from Eq. i|13a|) . The 
variational model is more accurate for Gaussian poten- 
tials. The main discrepancy we have found after a wide 
numerical exploration of the parameter space is a factor 
\[2 that is mainly due to the choice of the ansatz. 
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FIG. 3: [Color online] Effective potentials obtained for two 
different shapes of a(z): a step function and a Gaussian+step 
distribution. The potentials correspond to a = a cr for each 
form of the scattering length. The continuous line shows the 
effective potential for a = (the linear trap). The frame at 
the left correspond to the emission of a single soliton in the 
case of a Gaussian+step modulation of a (dotted curve). The 
shaded zone indicates the region of negative scattering length 
in the step model. 



N 




FIG. 4: [Color online] Same as Fig. EJfor the Gaussian+step 
distribution. The effective potentials correspond to different 
values of the penetration parameter D. 



IV. PARTIAL OUTCOUPLING AND 
MULTISOLITON EMISSION 



The distance D between the edge of the region in which 
the scattering length varies and the center of the trap 
plays an important role in the emission process as shown 
in Fig. 0] As it can be seen in the plots, there is a crit- 
ical overlapping between the a < region and the trap, 
for which the effective potential shows a central maxi- 
mum between two adjacent minima. In this case, part 
of the cloud is reflected backwards by the intermediate 
barrier and thus, it is not possible to extract the whole 
cloud (see inset of Fig. |3J. Thus, the emission will be of 
higher quality in traps of width below a certain critical 
threshold, which can be numerically calculated from the 
conditions aH/dz = and d 2 Tl/dz 2 > 0. 

A deeper numerical exploration based on Eq. re- 
veals interesting effects beyond those contained in the 
averaged Lagrangian description. An example is shown 
in Fig. [S] As a* becomes more negative we obtain the 
emission of an integer number of solitons. 

In this section and to approximate our analysis to re- 



alistic scenarios we have considered the scattering length 
to be of the form 



a(z) = a exp 



z-D 



(14) 



i.e., a supergaussian of finite width. 

In Fig. [S] we illustrate the dependence of the emission 
for different shapes of the function a(z). We begin with 
a profile close to the step form (m = 100) and then we 
relax the shape by diminishing the value of m until reach 
the Gaussian function. We have also observed that for 
finite a < regions, once the emitted solitons reach the 
opposite edge of the zone with negative scattering length, 
they are reflected backwards and thus remain trapped in 
the vicinity of the BEC reservoir [see Figs. OJd) and|Bl 
(a)]. This effect may have applications in the control of 
the emitted solitons once outcoupled and in the design 
of practical devices like laser tweezers for atoms [2(| . 

Another interesting consideration in order to construct 
experimentally wide regions of negative scattering length, 
is to use the superposition of several mutually incoher- 
ent laser beams as it shown in Fig. H3 In Fig. H3 (a) 
the region of negative scattering length is generated by a 
single outcoupling beam as in Fig[S]c). As a result of the 
finite width of the outcoupling beam, the atomic soliton 
rebounds at the boundary as is redirected backwards to 
the reservoir. In Fig. El(b) three equal Gaussian beams 
separated by its widths are used to generate a wider re- 
gion of negative scattering length. In the pictures, ver- 
tical axis is time form t = to t = lOOO^J 1 . The other 
parameters are the same as in Fig. |5]c). 

We have also observed in the numerical simulations 
that the soliton emission is more efficient if the scattering 
length varies abruptly in space. Moreover, for avoiding 
the return of the solitons once they have been emitted, 
it is desirable that the shape of the outcoupling beam be 
asymmetric. In the case in which laser beams are used 
to control the scattering length this could be achieved by 
partially overlapping parallel lasers. In order to generate 
a extended negative scattering length region, more lasers 
can be added or the outcoupling beam can be displaced 
following the motion of the emitted soliton. 



V. EFFECT OF TEMPORAL MODULATION OF 
THE OUTPUT 

Finally, we have also studied the temporal control of 
the emission process by simulating the dynamics of the 
solitons outcoupled with pulsed beams. We have run cal- 
culations for different temporal sizes of the beams. The 
results are shown in Fig. \7\ In these simulations we 
have used temporal profiles made with a Gaussian ramp 
connected to a flat top of variable size t p . The spatial 
profiles are step functions. As it can be seen in the pic- 
tures, depending on the duration of the pulse (i.e.: the 
size of the flat top), the quality of the outcoupling varies, 
yielding to a broadening in the profile of the solitons as 
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FIG. 5: [Color online] Emission of atomic solitons for different 
outcoupling laser beams (shaded region in the top plots) of 
supergaussian shapes with different values of the parameter 
m: a) m = 100, b) m = 8, c) and d) m = 2 (Gaussian 
distributions of different widths). The small frames in the 
top of each picture display the profile of the solitons at t = 
6001/J 1 . Vertical axis in bottom pictures is time from t = 
to t = 600z/J . Horizontal axis is 60 times the width L of 
the gaussain trap. In all cases the product Na is 8 times the 
value of Fig. 21 The rest of the parameters are the same as 
in Fig. 01 



the nonlinearity vanishes. In some cases the solitons can 
fuse after the emission. By varying the time that a is 
switched to a negative value, it is possible to obtain a 
periodic reconstruction of the individual solitons. The 
results are summarized in Fig[7|[c).In this simulation, the 
spatial profile of the outcoupling beam is a step function 
that is modulated in time with a Gaussian ramp con- 
nected at its maximum to a flat top of variable size that 
ends with a final Gaussian decay. The captions corre- 
spond to different sizes of the plateau t p and separation 
between applied pulses t s . In all cases time goes from 
t = to t = 500^2 : . a) corresponds to a continuous 
beam (i.e.: no modulation), in b) and c) we applied three 
pulses with t p = lOO^J 1 in both cases and t s = hv^ 1 and 
t s = Ibv^ 1 , respectively. The rest of the parameters are 
the same as in a). 



VI. CONCLUSIONS 

In summary, we have analyzed in detail the recently 
proposed mechanism for outcoupling coherent matter 
wave pulses from a Bose-Einstein Condensate. By us- 
ing this technique it could be possible to obtain a regular 
and controllable emission of atomic soliton bursts that 



7 





FIG. 6: [Color online] Emission obtained with the superposi- 
tion of parallel Gaussian beams. In (a) the region of negative 
scattering length is generated by a single outcoupling beam as 
in Fig|^|c) . In (b) three equal Gaussian beams separated by its 
widths are used to generate a wider region of negative scatter- 
ing length. Top plots display the output at t = lOOOz/J^ 1 . In 
the pictures, vertical axis is time form t = to t = lOOOi/J 1 . 
The other parameters are the same as in Fig. c). 
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FIG. 7: [Color online] Effect of a temporal modulation in the 
emission process. The spatial profile of the outcoupling beam 
is a step function that is modulated in time with a Gaus- 
sian ramp connected at its maximum to a flat top of variable 
size that ends with a final Gaussian decay. The captions cor- 
respond to different sizes of the plateau t p and separation 
between applied pulses t s . In all cases time goes from t = 
to t = 500UJ 1 . a) corresponds to a continuous beam (i.e.: 
no modulation), in b) and c) we applied three pulses with 
t p — 100^ 1 in both cases and t a = bv^ 1 and t a = 15i^J , 
respectively. The rest of the parameters are the same as in [5] 
a). 



are easily extracted by an adequate choice of the control 
parameters. 



The particular shape of the spatial dependence of the 
scattering length used for outcoupling the solitons does 
not affect essentially the emission process although it is 
of higher purity for sharp variations of the scattering 
length and for regions of negative scattering length of 
finite extent the soliton is reflected back when it reaches 
the boundaries. The temporal control of the solitons can 
be easily implemented with pulsed beams. 



We must stress that our results hold for different 
atomic species like 85 Rb and 133 Cs, with adequate pa- 
rameters. For the case of 133 Cs the scattering length can 
be controlled with high precision and can be done nega- 
tive and large. This could provide a system in which a 
very high number of solitons could be generated. 



Using the mechanism proposed in this paper a train 
of even several hundred of solitons could be coherently 
outcoupled from a condensate. As the techniques for co- 
herently feeding the remaining condensate progress our 
idea could provide an outcoupling mechanism for a con- 
tinuous atomic soliton laser. 
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